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Controlled generation and manipulation of photon states encoded in their spatial degrees of free- 
dom is a crucial ingredient in many quantum information tasks exploiting higher-than-two dimen- 
sional encoding. Here, we prove the impossibility to arbitrarily modify d-level state superpositions 
(qudits) for d > 2, encoded in the transverse modes of light, with optical components associated to 
the group of symplectic transforms (Gaussian operations). Surprisingly, we also provide an explicit 
construction of how non-Gaussian operations acting on mode subspaces do enable to overcome the 
limit d = 2. In addition, this set of operations realizes the full SU(3) algebra. 
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I. INTRODUCTION 

Most promising approaches for scalable quantum com- 
munication (QC) rely on the use of photons as the main 
carriers of information among remote nodes of quantum 
networks, where matter-based quantum memories are lo- 
cated [TJ |2j [3j HI [5] . Photons, besides being the natural 
candidate for QC due to their long decoherence time and 
the relative ease with which they can be manipulated, 
can actually encode multiple quantum bits of information 
(qubits) into various degrees of freedom. These include 
frequency, polarization, linear momentum and orbital an- 
gular momentum [B]. The possibility of simultaneously 
exploiting these degrees of freedom [Tj H] is becoming in- 
creasingly appealing for the faithful mapping of quantum 
states between light and matter [3] . A fundamental ques- 
tion then arises: What are the most general photon state 
manipulations allowed by benchmark optical components? 
It is of paramount importance to obtain a clear represen- 
tation of all such state mappings to further develop a 
truly multi-degree-of-freedom photon state engineering. 

Here, we address the problem of whether, by resort- 
ing to the symplectic group of optical transformations 
on spatial transverse-field modes, it is possible to per- 
form arbitrary manipulations on photon states encoded 
in large, but finite, d-dimcnsional superpositions of these 
modes (quefits). This is relevant for non-dichotomic 
QC protocols, which include those exploiting multimodc 
squeezing [9] and the orbital angular momentum (OAM) 
of light [lUl [TU HH [IHl IIH US! HZ]. For OAM, one 
of its main distinguishing features is the access to, in 
principle, an infinite-dimensional Hilbert space expanded 
by cylindrically-symmetric paraxial eigenmodes (e.g. the 
Laguerre-Gaussian basis) [TBI HH HOJ HI]- Spatial en- 
coding conveys several independent channels of informa- 
tion that could be very useful in quantum cryptographic 
schemes with larger alphabets [22] and security enhance- 
ment against eavesdropping 123] . Even for quantum com- 
putation applications, the high-dimensional aspect would 
enable to optimize certain computing arquitectures |24j . 

A necessary condition to perform arbitrary unitary op- 
erations on a pure quantum state \ip) = Ylj=i a j\j)' 



consisting of a d-dimensional superposition of orthogonal 
eigenmodes is to modify in a controlled way each of 
the complex coefficients Oj . In most of the experimental 
realizations oriented towards the use of spatial degrees of 
photons for high-dimensional encoding, phase holograms 
and reconfigurable spatial light modulators have been 
employed to approximately manipulate specific combi- 
nations of optical transverse modes [BJ. In practice, how- 
ever, these elements do not strictly preserve paraxiality 
but, rather, behave as non-unitary transformations, thus 
constituting a source of mode noise that eventually de- 
stroys the desired large, but finite, multidimensionality of 
the quantum states to be exploited. Our first main result 
shows that when these, or any combination of, optical el- 
ements belong to the group of symplectic transformations 
(which include Gaussian operations), it is impossible to 
arbitrarily modify single-photon qudit states for d > 2, 
via unitary operations generated by those transforms. 
Hence, a clear motivation emerges: Is it possible to find 
transformations on paraxial modes which allow one to re- 
ally overcome the limit d — 2? Our second main result 
provides a positive answer to this question; we present a 
set of non-Gaussian operations that truly enable us to ar- 
bitrarily manipulate (up to global phases) single-photon 
qutrit (d = 3) states. Furthermore, this set of operations 
constitutes a SU(3) algebra. 

The paper is organized as follows: Section [H] gives a 
brief summary of the formalism on symplectic groups and 
transformations in the optical phase space. In Section [TTT| 
we introduce and characterize the most general represen- 
tation of unitary (metaplectic) operators corresponding 
to all possible optical symplectic transformations that 
can be performed on transverse field modes. Section |IV| 
provides the first main result of our paper; we prove that, 
via the group of symplectic transforms acting on super- 
positions of paraxial modes, it is impossible to imple- 
ment operations that change arbitrary qudit states onto 
any other qudit for d > 2. In Section [V] we extend our 
analysis to non-Gaussian operations on these modes and 
present new routes towards the aim of truly manipulat- 
ing arbitrary single-photon qudits. Section |VI| concludes 
the paper with a discussion of alternative approaches to 
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implement controlled gates on single-photons using more 
than one of their degrees of freedom. A simple optical 
scheme for a CNOT gate exploiting OAM and polariza- 
tion, is proposed. 



II. SYMPLECTIC GROUP FORMALISM 

To put in context the class of optical transforma- 
tions referred to above, it is necessary to start by in- 
troducing the symplectic formalism that will be used 
extensively throughout the paper. We first recall that 
the dynamics of classical and quantum Hamiltonian sys- 
tems has an underlying symplectic structure. Symplec- 
tic methods have been applied in the theory of elemen- 
tary particles, condensed matter, accelerator and plasma 
physics, oceanographic and atmospheric sciences and in 
optics [551 1551 [57] ■ Fundamental to all of them is the 
phase space picture. Any classical system with n de- 
grees of freedom is described by a set of pairs qj,Pj 
(j = 1,2, ... ,n) of mutually conjugate canonical vari- 
ables. In the quantum domain one can associate to 
these variables the irreducible set of canonical Hermi- 
tian operators qj,Pj- The basic kinematic structure is 
provided by Poisson brackets in the former case and by 
the Heisenberg commutation relations in the latter. By 
assembling the canonical variables and operators into 2n- 
component vectors £ = (gi , q 2 , . . . , q n , pi , P2 , ■ ■ ■ , Pn) an d 
i = {qi,q~2,---,qn,Pl,P2,---,Pn), the Poisson brackets 
and the Heisenberg commutation relations can be cast, 
respectively, as {£a,£,a} = n 0)/3 and [i a ,£p] = i^a,/3, 
(a,/3= 1,2,..., 2n), where 
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is the 2n-dimcnsional symplectic metric matrix. Of par- 
ticular relevance are the real linear canonical transfor- 
mations among quantum (classical) canonical quanti- 
ties [55]. They preserve the Heisenberg (Poisson) rela- 
tions and are represented by symplectic matrices S : £ — > 

£ = S£, obeying the condition SQ,S = £1. The set of all 
such 2n-dimensional real matrices forms the (2n 2 + n)- 
parameter non-compact symplectic group Sp(2n,R). 

The power of symplectic formalism becomes apparent 
in the following general setting. Let Ti denote the Hilbcrt 
space of n-mode states p on which £ Q act. Given that for 
any S € Sp(2n, R) the Hcrmiticity properties and com- 
mutation relations of the £ Q are conserved, and since £ a 
act irreducibly on Ti., it follows from Stone- von Neumann 
theorem that one can define unitary operators U(S) on 
Ti, implementing SQS T — Q, such that (28] 



(2) 



Clearly, all possible transformations S are mapped, up to 
a sign ambiguity, onto U(S). Hence, all density operators 
pen transform under U(S) as p' = U{S)pW{S). On 
the level of Wigner functions W(£), this transformation 
acquires a strikingly simple form: W'(£) = t¥(5' _1 £). 



III. METAPLECTIC OPERATIONS ON 
PARAXIAL MODES 

Within classical and quantum optics, the symplectic 
formalism has been extensively used both in studying 
mode-mapping properties of lossless first-order (paraxial 
or ABCD) systems [HJ HOI EE [301 131] and in char- 
acterization of continuous- variable entanglement [32 . 
An important class of symplectic transforms is that of 
two-mode systems represented by the symplectic group 
Sp(4, R). For instance, bipartite Gaussian operations, 
which preserve the Gaussian character of the Wigner 
functions, belong to Sp(4,R). Let us make explicit the 
form of all possible U(S) when S S Sp(4, R). They give 
rise to the unitary metaplectic representation of Sp(4,R) 
acting on H. All these unitary operations U(S) are gen- 
erated by ten Hermitian operators ^7", quadratic in £, 
that can be split into two sets [33]: passive and active 
generators. The passive set encompasses the maximal 
compact subgroup U(2): 
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Here, cij — (qj + if>j)/\/2, (respectively a|), j — x,y, 
are the two annihilation (creation) operators for all or- 
thogonal transverse modes. The passive operators ^ 
have the form of the well-known Stokes operators. They 
obey the usual commutation relations [Ci,Cj] = ietjkCk 
k — x,y, z), with C a being the only commuting ele- 
ment in U(2). 

The active set is responsible for the noncompactness 
of Sp(4,R): 
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satisfying the following commutators: 

[Ci,K.j] = ie ijk ICk, [£i,Mj] = ieijkMk, 
[t ,Kj ± iMj] = ± iMj) , [KuMj] = i , 
[/Ci,/Cj] = [Mi,Mj] = -ieijkCk- 



Below we elucidate the action of each of the ten genera- 
tors on a spatial mode carrying OAM. 

As any arbitrary sequence of symplectic transforma- 
tions S m is again another symplectic transformation S = 
Y[ m S m , one concludes that the most general metaplec- 
tic operator U (S) corresponding to S is represented by a 
single exponential of i times real linear combinations of 
any of the above generators: 



U(S) = exp 



(-is ■ j) 



(5) 



with J~ £ {£, ft, Ai}, and s a ten-parameter vector. 

When applied to photon number states, the passive 
(active) generators have a well-known interpretation: 
they conserve (do not conserve) photon number. Un- 
like active generators, which require nonlinear photon in- 
teractions, passive generators can be implemented with 
linear optical components: beam splitters and phase 
shifters. Now, despite the exact isomorphism between 
symplectic transformations on photon number states and 
spatial modes, they have quite distinct physical implica- 
tions. To gain insight on how the metaplectic opera- 
tor (|5| affects spatial modes, we resort to the Wigner 
representation in conjunction with the Stone-von Neu- 
mann theorem ^ . A revealing example is the following. 
Consider a Laguerre-Gaussian mode [B] LG^ jP , where the 
indices I — 0, ±1,±2, ... and p — 0,1,2,... stand for 
the topological charge and the number of nonaxial radial 
nodes. Let We tP (£) be the associated Wigner function, 
which, in the general case, is non-Gaussian jTSj 
For each of the ten generators ^ and (EJ one can eas- 
ily obtain the corresponding symplectic matrix S acting 
on the input Wf iP (£) to yield exact analytical expres- 
sions for the output Wigner functions (via its covariant 
property). The resulting position distribution i| '(q) is 
computed via the marginal jg(q) = / (PpW^S^i). 

Figure [l] depicts /jf^ p=0 (q) under the action of each 
passive and active generator. It can be seen that they 
produce fundamentally different mode-mapping geome- 
tries. Passive generators describe rotations on the or- 
bital Poincare sphere [T5J [201 131] • They preserve the 
order N = \£\ + 2p of any of the modes lying on the 
sphere. Generator C yields the mode-order, C z repre- 
sents real spatial rotations on the transverse x — y plane 
containing the modes and is proportional to the compo- 
nent of the OAM operator along the propagation direc- 
tion |20j . with LG^p being the eigenmodes. C x and C v 
represent simultaneous rotations in the four-dimensional 
phase-space: C x produces rotations in the x — p x and 
y — p y planes by equal and opposite amounts, whereas 
C y gives rise to rotations in the x — p y and y —p x planes 
by equal amounts. The cigcnstates of C x are the Hermite- 
Gaussian HG„^ jnH modes, where n x ,n y are nonnegative 
integer indices, and their mode-order is N = n x + n y . 
Both Laguerre- and Hermite-Gaussian bases are unitar- 
ily related: LG^p transforms into HG„.. m „ via e~ 1 ^ / 2 ^ Cy 
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FIG. 1: Action of passive (a), and active (b) generators on 
a Laguerre-Gaussian mode LG^=i jP= o. Profiles in (a) are 
mapped onto the first-order orbital Poincare sphere. Profiles 
in (b) are given for increasing squeezing parameter s. 



We note in passing that the interferometric scheme pro- 
posed in Ref. [3S] to measure the OAM spectrum (i.e. 
the i index-decomposition) of a light beam, which only 
involves C z , could be generalized to determine the com- 
plete Hermite-Gaussian spectrum by replacing C z with 
C and C x |36j . In contrast with passive generators, the 
active ones scale {squeeze) the spatial modes and change 
the order N , giving rise to infinite mode-superpositions. 
Of these, only ft z and M y suffice to describe, jointly 
with the set the general metaplectic operator ^ 
by recourse to the following passive-active-passive de- 
composition U = e -i^C e -i{v y M y + v ^) e -ir 1 C i sti n rc _ 
quiring ten parameters. The symplectic matrices asso- 
ciated with both passive and active generators can be 
implemented with a small arrange of (< 10) spherical 
and/or cylindrical lenses solely controlled by variations 
of the focal lengths and/or rotations along the system 
axis [nni SZ1 |3S]; that is, with simple linear optical com- 
ponents. 



IV. METAPLECTIC OPERATIONS ON 
SINGLE-PHOTON QUDITS 

Having characterized the complete set of unitaries gen- 
erated by passive ^ and active (|4| operators that can 
be performed on paraxial field modes, we proceed with 
all possible actions of the metaplectic operator ^ on 
arbitrary single- photon qudit states. The most gen- 
eral (paraxial) single-photon pure state can be described 
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as [20] 

p OO 

|V) = Yl / ^C-,n«,n B H6U,n B M|vac) . (6) 
cr,n K ,n y ^ 

Here, n<o n (w) denotes the bosonic creation opera- 
tor of a Hermite-Gaussian mode, linear polarization a, 
and frequency u). The commutation relations read as 

[ba,n x ,n v M ,^', < ,„; (^')] = "W^n^n^C^ - fc/). 

The normalized complex coefficients C CT .„ x .„ y (cj) can be 
interpreted as the probability amplitudes for finding a 
photon in the state b\ n ^ n ^ (lu) |vac) = |o")0|tIx, %)<8|w). 
Let us concentrate on the spatial part of Eq. ^ and as- 
sume that it consists of a finite superposition of d or- 
thogonal HG modes \n xi n y ) = (a*)™* (at) n * |0, 0) (|0,0) 
is the fundamental Gaussian mode), in the normalized 
quefit form 

W) = c n*,n y \n x ,n y ) . (7) 

n x ,n y 

Any qudit requires, at least, 2d independent real param- 
eters, albeit normalization and invariance of Q under a 
global phase reduces this number to 2(d— 1). A necessary 
condition to fully manipulate a single-photon state ^ is 
to arbitrarily modify the d complex coefficients c n ^ ny 
(e.g. it should be possible to set all coefficients c„^„ y 
equal to zero, except for one of them), leaving invariant 
the d-dimensional subspace Hd expanded by {\n x ,n y }}. 
In other words, one must discard all transformations 
on ^ giving rise to modes not belonging to Hd- 

Let us analyze the most important restrictions imposed 
by unitary operations acting on ([T]) and generated by 
the group of transformations S <E Sp(4, R). First, no- 
tice that since the general metaplectic operator ^ in- 
volves ten generators, ci-dimensional superpositions with 
d > 6 cannot be arbitrarily transformed within Sp(4, R). 
This fact, of course, does not preclude the possibility 
to manipulate qucfits with d < 6. A key observation is 
the recognition that finite dimensional representations of 
Sp(4,R) are necessarily nonunitary, owing to the non- 
compactness of Sp(4,R). That is, the noncompact part 
of Sp(4,R), represented by the active generators Q, is to 
be excluded from the set of symplectic transformations 
in order to maintain the subspace Tid finite. Otherwise, 
the qudit ^ would become an infinite superposition of 
all HG modes under the general action of More ex- 
plicitly, let Vd be a projector in Hd, so that it fulfills 
"PdW) — \4>)- Notice that if the metaplectic operator ^ 
must keep Q in Hd, then PdU\ip) — U\ip), which im- 
plies that U~ lr PdU — Vd- This last condition should 
hold for any choice of the parameters s in (|5j, and so it 
follows that the commutator [s • ij, Vd] — 0. However, 
this vanishing commutator is incompatible with the pres- 
ence of the active generators Q. One is therefore lim- 
ited to the compact subgroup of Sp(4,R), i.e. to the set 
of passive generators ^ , to perform unitary operations 



on ([7]) leaving the subspace Hd invariant. This means 
that the metaplectic operator (|5| reduces to the passive 
one Ujc — e - l i s " c "+s :l: c x j rSyCy+s,c,) ^ which, now, only 
contains the four free parameters {s , s x , s y , s z }. Conse- 
quently, one must confine to three-dimensional subspaces 
Hd, with state ^ being a qutrit. In principle, the in- 
tervening modes in ^ can have different order N. How- 
ever, since Uc preserves N, this implies that it would 
then be impossible to attain the mode transformation 
\n x ,n y ) — ► \n' x ,n' y ) when Ti x ~\- TLy ^ n' x +n'. There is still 
the possibility that the three modes could have the same 
order. In this particular case, taking into account that 
C \ip) = N/2\tp) and = 0, one can express the 

action of U c as U c \^) = e -«o JV / 2 e -*( s *A*+ s vA'+^*)|^). 
Up to a global phase, there are only three free parame- 
ters {s x , s y , s z } to carry out the general transformations 
on the state ([7]), which are insufficient even for qutrits 
(as they involve, at least, four free real parameters). We 
have thus proven the following result: 

Proposition.- It is impossible to arbitrarily modify the 
c?-dimcnsional mode-superpositions of single- photon pure 
qudit states ^ for d > 2, via unitary operations U(S) 
generated by symplectic transforms S E Sp(4,R). 

Several comments are in order. Our proposition is ex- 
pected to also hold for mixed states. Arbitrary opera- 
tions on qubits (d = 2) are not prohibited within the 
subgroup U(2) of Sp(4,R), as it should be [20! ■ Operators 
Uc acting on higher-than-two mode superpositions ^ re- 
strict the possible values of coefficients c„ x <n ; the higher 
the value of d the larger the number of constraints on 
Cn m ,n y - hi fact, it is easy to show that, within Sp(4,R), 
the most general finite superpositions ^ reduce to the 
well-known Dicke or spin coherent states |39j (which de- 
pend only on two real parameters, leaving aside global 
phases) . Although we have identified spherical and cylin- 
drical lenses as the main optical elements of Sp(4,R), the 
above proposition also affects phase holograms belong- 
ing to the metaplectic representation of Sp(4,R). Not all 
unitary and paraxial transformations are in such repre- 
sentation. It is an open question to determine the entire 
class of unitary operators outside the metaplectic rep- 
resentation of Sp(4,R) that leave invariant the bases of 
paraxial modes; but it would definitely fall into the cat- 
egory of non-Gaussian operations. In the next Section 
we partially clarify this question; we find examples of 
non-Gaussian operations which preserve paraxiallity and 
allow us to overcome the limit established by the preced- 
ing proposition. Moreover, since Laguerre- and Hermite- 
Gaussian bases are unitarily connected, our proposition 
also establishes the impossibility to achieve arbitrary 
quefit gates on multi-dimensional superpositions of modes 
bearing OAM. Given that the concept of OAM of light 
is only strictly meaningful within the paraxial approxi- 
mation [5] [20] 1 encoding photons in non-paraxial modes 
would generally couple polarization and OAM making 
most QC tasks in such scenarios extremely difficult. 



V. NON-GAUSSIAN OPERATIONS 

In contrast with Gaussian operations, non-Gaussian 
operations remain to be fully explored. It has been rec- 
ognized that non-Gaussian operations could represent an 
advantage to perform some quantum information tasks. 
Non-Gaussian operations on continuous variables allow 
the access to results beyond no-go statements concerning 
Gaussian operations. For instance, quantum speed up is 
impossible for harmonic oscillators by Gaussian opera- 
tions with Gaussian inputs |40j . Distillation of Gaussian 
bipartite entanglement is also impossible by performing 
only Gaussian local operations and classical communi- 
cation based on homodyne detection and requires non- 
Gaussian operations |41j . It has been experimentally 
demonstrated that entanglement between Gaussian en- 
tangled states can be increased by conditional subtrac- 
tion of single photons from the Gaussian beams [42] . For 
universal quantum computation with continuous vari- 
able cluster states [12] it is necessary, at least, one non- 
Gaussian projective measurement. 

In the present scenario, our preceding proposition im- 
poses a restrictive limit for Gaussian operations on arbi- 
trary superpositions of spatial mode states. However, in 
this section, we report a new class of non-Gaussian oper- 
ations which enable us to fully manipulate superpositions 
of three-level states (qutrits) , beyond the restrictions im- 
posed by our above no-go proposition. Furthermore, this 
class of operations forms a complete set of single qutrit 
gates fulfilling a SU(3) algebra. 

Consider the eight following generators acting on the 
Hilbert space of Hermite-Gaussian modes: 

% = -r (ola y - a+o x ) , 
% = o - a\a y ) , 

^4 = 2 ("i + ^ - a x a x a x - a x a x a x 

- a' x a y a y - a x a y a y ) , 

% = --(a x - a x - a) x a x a x + a x a x a x 

- a x a y a y + a x a y a y ) , 

% = 2 (°i + ®v ~ tiyQyG'y ~ o-ld y a y 

- a\ala x - a y a\a x ) , 

% = - - (a+ - a y - a y a y a y + a y a y a y 

- a y a\a x + a y ala x ) , 

% = [-2 + 3 +0+0,,)] . (8) 

Quite remarkably, these generators, within the sub- 
space expanded by the Hermite Gaussian modes Ht = 
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|1,0) |0,1) 

FIG. 2: Mode conversion via the unitary operations Ur gen- 
erated by the SU(3) group 

{|0,0),|1,0),|0,1)}, obey the SU(3) algebra [f a ,%] = 
ifabcT~ c (a, b, c = 1, 2, . . . , 8), where the only nonvanish- 
ing (up to permutations) structure constants f a bc are 

given by / 123 = 1, /l47 = /l65 = /246 = /257 = 

/345 = /376 = 1/2, and / 458 = f e78 = V3/2. Un- 
like in the metaplectic representation of Sp(4,R), there 
are now four generators in Q involving cubic terms in 
a x ,a x ,a y ,a y . Their structure resembles that of fully- 
quantum interaction Hamiltonians found in nonlinear 
parametric processes. Notice that the triad of genera- 
tors Ti = {7i,72,73} corresponds to the SU(2) group 
in Eq. ([3]). The remaining two SU(2) groups are formed 
by the triads: T 2 = {%,%, (f 3 + V3T 8 )/2}, and T 3 = 
{%,7~7, (—73 + v / 37s)/2}. Unitary operators Ur 1 gener- 
ated by the first triad give rise to superpositions between 
the two modes |1,0) and |0, 1), leaving invariant the fun- 
damental Gaussian mode |0,0). Unitaries Ur 2 and Ur 3 , 
generated by the second and third triad, produce super- 
positions between the two modes |0, 0) and |1, 0) (leaving 
invariant |0,1)), or the modes |0, 0) and |0, 1) (leaving 
invariant |1,0)), respectively. In contrast with Ur 1 , the 
action of Ur 2 an d Ur s on gives rise to a new fea- 
ture: non-conservation of the mode order. However, this 
non-conservation is fundamentally different to the one en- 
countered in the noncompact representation of Sp(4,R), 
since it preserves the subspace Ti-r . Figure [2] summa- 
rizes the action of the three unitaries Uj< on the subspace 
Ht- It is worth mentioning that all these non-Gaussian 
operations can be implemented with passive optical el- 
ements having higher-than-first-order aberrations (non- 
quadratic refractive surfaces) |44j . An open and inter- 
esting problem would be the extension of the Stone-von 
Neumann theorem ^ to the case of our qubic genera- 
tors ^ . This would enable one to find the explicit form 
of the symplectic transform and thus the construction of 
the associated optical system. 

Complete manipulation of the qutrit = co,o|0, 0) + 
ci,o|l,0) + cq,i|0, 1) is now possible using the SU(3) 
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FIG. 3: (Color online) Scheme of a single-photon polarization- 
CNOT gate. According to the polarization of the in- 
put Laguerre- Gaussian mode the photon experiences a sign 
change (£ -> —£, where \£\ = 1) in its OAM. 



group ^ , although to produce a general qutrit (up to a 
global phase) it suffices to perform the following sequence 
of operations: starting, for example, with an input fun- 
damental Gaussian mode |0, 0) subjected to the unitary 
operator Ur 2 , one can obtain 

#r a |0,0)=coB~|0,0) + e*sin~|l > 0>. (9) 

Then, taking into account the closed SU(2) algebra 
obeyed by Urn it follows that 



C/r 1 |l,0) = cos-|l ! 0) + e^ sin-|0,l). (10) 



With this specific operation structure, using (JO]) and ((10), 
we can construct a general normalized qutrit state (up to 
a global phase) : 



[/ ri C/r 2 |0,0) = cos-|0,0) + e^sin-cos-|l,0) 

+ e l{ip+v>) sm e -sm 6 -\0,l) . (11) 

Since the four parameters d,9',tp,(p' can be varied in- 
dependently during the process, not all generators in the 
set (|8| are actually needed to produce any qutrit encoded 
only in paraxial spatial modes. Notice that our results 
can also be extended to other physical scenarios, due to 
the isomorphism of the formalism, although in our case 
an additional motivation is provided by the simplicity of 
their experimental implementation. 



VI. DISCUSSION AND CONCLUSIONS 

In spite of the stringent limits raised by our results 
of Section |IV[ we have clearly shown how non-Gaussian 
operations can circumvent most of those difficulties. It 
is also worth emphasizing that there exist other alterna- 
tive approaches, exploiting the spatial encoding of light, 



to fully manipulate higher-than-two-dimensional Hilbert 
spaces for various quantum information tasks. These ap- 
proaches rely on the use of several degrees of freedom 
of light, albeit they cannot attain very large subspace- 
dimensionalities. To illustrate, consider that we perform 
transformation (10) and, analogously to transformation 



<M, we wish to produce a general (up to a global phase) 
single-photon qutrit state. Instead of ([9|, we can use 
another degree of freedom of the same photon (e.g. po- 
larization). To do so, we need both a complete set of 
single-photon qubit gates in each of the degrees of free- 
dom and a conditional gate between the two involved de- 
grees of freedom. For instance, an efficient linear single- 
photon CNOT gate in which photon polarization acts 
as the control-qubit on the other photon degree of free- 
dom, OAM, which plays the role of the target-qubit, is 
feasible with current technology. There are several possi- 
ble routes, such as with space-variant optical axis phase 
plates made of nematic liquid-crystals [15] H3], or using 
a Mach-Zender configuration [47]. A conceptually sim- 
ple scheme is depicted in Fig. [3] This gate includes one 
polarizing beam splitter (PBS) and two pairs of cylin- 
drical lenses (CL) whose bases subtend a 45° angle with 
the plane of the interferometer. This interferometer re- 
sembles previous Sagnac interferometers used for mea- 
suring the spatial Wigncr function [15] and for other 
single-photon quantum gate demonstrations employing 
polarization and continuous variables [49]. There, the 
inner arm of the interferometer contained a Dove prism. 
Here, the presence of cylindrical lenses constitutes the 
key feature to exploit photon OAM. According to the 
input photon polarization state (horizontal or vertical), 
the input photon views the cylindrical lenses CL1 with a 
different orientation and experiences a mode transforma- 
tion (LG — ► HG) depending on the particular value of 
I. After exiting through the PBS, the second cylindrical 
lenses CL2 yields the mode transformation HG — > LG., 
which completes the action of the polarization-OAM- 
CNOT single- photon gate. 

An even more fascinating scenario is the transfer 
of photons carrying OAM onto Bose-Einstein conden- 
sates [50] and their storage in electromagnetically in- 
duced transparency media [51] . In this respect, it would 
be very interesting to explore the possibility of mapping 
the correlations of photons entangled in OAM [10] [12] 
IT51 ITU ITrJl HTfl |2"T] on quantum holograms, allowing for 
the reconstruction of nonclassical states of light from a 
matter-based quantum memory. 

In conclusion, we have shown that if single-photon pure 
qucfit states are prepared in a d-dimensional superpo- 
sition of spatial modes, it is impossible to arbitrarily 
change such mode-superpositions for d > 2, by solely 
resorting to unitary operations generated by symplectic 
transforms of the group Sp(4,R). Our results provide a 
complete characterization of linear canonical transforma- 
tions on transverse optical modes and pose a considerable 
challenge on quantum communication protocols exploit- 
ing multidimensional spatial encoding: one cannot have 
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a full access and control of large but finite-dimensional 
Hilbcrt spaces expanded by these modes. Implementa- 
tion of a new class of paraxial non-Gaussian transforma- 
tions to fully encode any arbitrary single-photon qudit 
state is required. We have provided an explicit construc- 
tion of this new class of operations. Moreover, using the 
spatial encoding in combination with other degrees of 
freedom, one can overcome this problem though at the 



price of scalability. In this conditional gate be- 

tween the involved degrees of freedom is needed. 
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